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ТЕОРИЯ РАСКЛИНИВАЮЩЕГО ДАВЛЕНИЯ 
В КОНЦЕНТРИРОВАННОЙ ДИСПЕРСИИ ЗАРЯЖЕННЫХ 
СФЕРИЧЕСКИХ ЧАСТИЦ 


1. ЯЧЕЕЧНАЯ МОДЕЛЬ ДЛЯ РЕШЕНИЯ УРАВНЕНИЯ ПУАССОНА — БОЛЬЦМАНА 


Жарких Н. И., Шилов В. Н. 


В работе предложен ячеечный метод вычисления уравнения состояния 
дисперсии (зависимости расклинивающего давления от объемной доли ча- 
стиц), стабилизированной отталкиванием перекрывающихся двойных слоев. 
Обсуждается соотношение ячеечного метода с другими способами вычисле- 
ния расклинивающего давления. - 


Постановка задачи. Традиционная постановка задач в теории устойчивости коллои. 
дов [1—4] такова: две коллоидные частицы сталкиваются в процесс броуновского двн- 
жения. Необходимо вычислить силу взаимодействия между этими частицами и вероят- 
ность образования агрегата при столкновении. Эти задачи возникают при анализе агре- 
гативной устойчивости разбавленных дисперсий. Другой тип задач связан с вычислением 
расклинивающего. давления и анализом устойчивости в концентрированных дисперсных 
системах. Задачи такого типа возникают при изучении периодических коллоидных струк- 
тур [5], динамических мембран [6]. 

В концентрированной дисперсии сила взаимодействия между двумя частицами силь- 
но зависит от расположения соседних частиц. Это означает, что нельзя ввести некий 
универсальный потенциал взаимодействия двух частиц и полную потенциальную энергию 
дисперсии представить в виде суммы двухчастичных потенциалов по всевозможным 
парам частиц (подобно тому, как выражается потенциальная энергия газа молекул). 
Поэтому наиболее важной характеристикой устойчивости концентрированной дисперсии 
является уравнение состояний — функция, связывающая расклинивающее давление в си- 
стеме коллоидных частиц с их объемной долей. В отличие от потенциала или силы взаи- 
модействия двух частиц уравнение состояния (иначе его можно назвать изотермой) — 
статистическая характеристика системы, заведомо содержащая не две, а очень большое 
‚количество частиц. 

В данной работе мы рассмотрим систему одинаковых заряженных сферических ча- 
стиц радиуса (), логруженных в водный раствор 1—1-валентного электролита с концен- 
трацией Со. Пренебрежем силами притяжения между частицами (силами Ван-дер-Вааль- 
са). Тогда частицы будут только отталкиваться за счет перекрытия их диффузных атмо- 
сфер. Если дисперсию заключить в некоторый заданный объем, частицы в ней будут 
стремиться разойтись на максимальные расстояния. В результате частицы будут рас- 
полагаться в узлах кубической гранецентрированной решетки (или эквивалентной ей 
в энергетическом отношении гексагональной решетки) и образуется периодическая кол- 
лоидная структура (ПКС) с ограниченным объемом [5]. Поскольку между частицами 
действуют только силы отталкивания, то с ростом объемной доли частиц (обозначим 
ее ©) расклинивающее давление будет расти и полученная ПКС будет термодинамически 
устойчива. Учет теплового движения частиц ведет к разрыхлению ПКС: если заданный 
объем достаточно велик (а расклинивающее давление достаточно мало), то понижение 
свободной знергии при образовании ПКС может быть так мало, что тепловые флуктуа- 
ции смогут существенно нарушить дальний порядок в решетке («плавление» ПКО). 
В данной работе мы, однако, пренебрежем влиянием теплового движения на свойства 
концентрированной дистерсии. 

Общие формулы для расклинивающего давления. Поскольку мы хотим вычислить 
давление, связанное только с электростатическим взаимодействием частиц, то необходи- 
мо найти только ту часть свободной энергии дисперсии, которая обусловлена зарядом 
частиц. Это просто сделать, рассматривая процесс зарядки первоначально нейтральной 
дисперсии от внешнего источника тока (рис. 1). Все частицы, находящиеся внутри сфе- 
рической оболочки, соединены с одним полюсом источника. Если при замыканий ключа 1 
(ключ 2 разомкнут) в цепи переносится заряд 49, то изменение свободной энергии со- 


<тавит 
ЧЕ=ФаО 


где ф — потенциал частиц. (Объем системы в процессе заряжения не меняется.) Полная 
хвободная энергия, связанная с зарядом, равна 


[#1 
Е= | 940 (1) 
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где О, — полный заряд, сообщенный системе. : 
Чтобы вычислить давление р, необходимо разомкнуть ключ /[, замкнуть ключ 2, под- 
ключающий схему, обеспечивающую характерную для данной системы связь заряда с 
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потенциалом. Назначение этой схемы в нащем мысленном эксперименте — смоделировать 
двойной слой, находящийся при постоянном потенциале (тогда схема 3 должна пред- 
ставлять собой источник напряжения) при постоянном заряде (тогда схема 3 представ- 
ляет собой просто разрыв в цепи) или какой-либо другой механизм формирования двой- 
ного слоя (например, соответствующий заданной константе диссоциации ионогенных 
грунн частицы; в таком случае схема 3 должна быть сложным нелинейным элементом). 
Воспользуемся формулой для дифференциала свободной энергии [7] 


которую можно переписать в виде 


дЕ 99 
р у Ри | (2) 


Здесь У — объем дисперсии, Е — свободная энергия, вычисленная по формуле (1). Фор- 
мула (2) — общее термодинамическое соотношение, позволяющее рассчитать давление 
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Рис. 1 Рис. 2 


Рис. 1. Эквивалентная схема концентрированной дисперсии заряжен- 
ных частиц, используемая при выводе формул для свободной энергии 
и расклинивающего давления 


Рис. 2. Диаграмма связи заряда дисперсии с ее потенциалом, исполь- 
зуемая при анализе термодинамической устойчивости ДС (на рис. 1,2 
пояснения в тексте) 


при любом способе распределения зарядов в системе (в частности, наличие второго чле- 
на в (2) позволяет вычислить р, если при изменении объема дисперсии меняется ее 
заряд). 

Каждый конкретный механизм формирования двойного слоя (ДС) приводит к опре- 
деленной зависимости заряда от потенциала: 


9=9(, и, .) 9 


для нахождения которой необходимо решить уравнение Пуассона — Больцмана с опре- 
деленными граничными условиями. 

Общие термодинамические соображения позволяют наложить некоторое ограничение 
на вид функции (3). Будем откладывать по осям координат рис. 2 заряд и потенциал, 
которые изменяются с изменением неких других характеристик системы. (например, 
объема, концентрации и состава электролита}. Тогда любая кривая вида (3), соответ- 
ствующая реальному механизму формирования ДС, должна лежать в пределах заштри: 
хованных квадрантов (в которых емкость се=@0/4ф>>0). Если дифференциальная ем- 
кость 49/4ф отрицательна, то система термодинамически неустойчива [7, $ 21] (элек- 
трическая емкость, как и теплоемкость,— вторая производная свободной энергии, кото- 
рая в тепловом равновесии должна быть положительна}. Края заштрихованной области 
на рис. 2 соответствуют постоянству потенциала (прямая 1) или постоянству заряда ` 
(прямая 2). Поэтому эти режимы формирования ДС можно считать предельными слу- 
чаями для всех возможных видов уравнения (3). Далее мы будем изучать этн предель- 
ные случаи, 

В случае 9О=сопз{ 


р=— а? 
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в случае ф==соп$ 


дЕ 00 
Е — 5 
р зу +5 (5) 


Выражения (1), (3)—(5) позволяют вычислить уравнение состояния дисперсии. 
Явный вид уравнения состояния зависит от модели ДС. В дальнейшем будем считать, 
что ДС имеет чисто диффузное строение. В таком случае для получения зависимости 
вида (3) необходимо решать уравнение Пуассона — Больцмана: 


} р ь 
ых ое (ехр (— $) — ехр(Ф)) (8) 


где Со нормальная концентрация электролита, =-- диэлектрическая проницаемость 
растворителя, Ф — потенциал ДС (символом 
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будет обозначать поле электростатического потенциала в отличие от величины р — зна- 
чения этого потенциала на поверхности частицы). 
Краевые условия для уравнения (6) имеют вид 


9Ф 
$ 
аФ 
в(®. 3%.) =° ® 
5; 


{$ — поверхность объема, ограничивающего дисперсию, $; — поверхность {-й частицы). 
Условие (7) обеспечивает электронейтральность дисперсии в целом, условие (8} — связь. 
между зарядом и потенциалом поверхности, характерную для данного механизма обра- 
зования ДС. В качестве явного вида функции [1 будем использовать либо 


Ф 15, = ф = с0п34 (9) 
либо 
0Ф 
— и (10). 
дп 5 248 


Ячеечный метод решения уравнения Пуассона — Больцмана. Для решения краевой 
задачи (6), (7), (9), (10) можно использовать различные методы. Например,. можно. 
разложить потенциал в ряд Фурье — Бесселя по сферическим гармоникам и вычислять. 
несколько первых членов этого ряда. Этот метод даже в дебаевском приближении при- 
водит к очень громоздким выкладкам; кроме того, необходимость усечения ряда при-- 
водит к ограничению области применимости такого подхода большими значениями хй 
(х — обратная дебаевская длина, # — зазор между частицами; практически хй>1). 
В данной работе для решения уравнения (6) предлагается ячеечный метод, позволяю- 
щий значительно упростить. вычисления и получить выражение, аппроксимирующее- 
истинное уравнение состояния (по крайней мере качественно} при всех значениях объем- 
ной доли и сильном перекрытии двойных слоев соседних частиц. Сильное перекрытие: 
в данном случае означает, что могут перекрываться дебаевские атмосферы не только бли- 
жайших, но и удаленных соседей. В случае тонкого ДС не сильное перекрытие в ука- 
занном смысле невозможно; тут ячеечная модель неприменима и следует пользоваться’ 
теорией ДЛФО 11, 2]. 

В применении к данной задаче идея ячеечного метода состоит в том, что сложную 
для вычислений систему краевых условий (7), (9) или (7), (10) можно заменить такими: 
же по виду краевыми условиями на поверхности каждой элементарной ячейки ПКС (т.е. 
под поверхностью $ в (7) следует понимать поверхность элементарной ячейки). Тогда‘ 
отдельные ячейки становятся независимыми и найдя потенциал внутри одной ячейки, 
можно продолжить его периодически (с периодом решетки ПКС) на весь объем диспер- 
сии. Физически такая замена краевых условий оправдана тем, что в ИКС все ячейки. 
одинаковы и условие электронейтральности выполняется ие только для дисперсии в це- 
лом, но и для каждой ячейки отдельно. Поэтому для ПКС такая модификация краевых 
условий является точной. 

Далее заменим приближенно точный многогранник элементарной ячейки ‘сферой с 


объемом, равным объему ячейки (радиус этой сферы равен 6==аа”\? ). Тогда краевые: 
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‘условия примут вид: 


оф 
ды =0 
5 | (10) 
Ф| ._«=Ф (12) 
или 
| __ 9 ы 
дг а  эАТа (18) 


Решение краевой задачи (6), (11)—(13) уже не представляет большого труда. В по- 
следующих работах мы найдем это решение для малых (приближение Дебая) и умерен- 
ных (приближение Гронвелла — Ла-Мера — Сандведа) поверхностных потенциалов и 
обсудим полученные уравнения состояния. 

В заключение отметим, что возможен иной подход к вычислению уравнения состоя- 
ния — прямой расчет конфигурационного интеграла системы при помощи Монте-Карло 
или какого-либо приближенного метода [8]. При таких расчетах надо исходить из неко- 
го явного выражения для потенциала взаимодействия двух частиц. Как отмечалось во 
введении, при сильном перекрытии ДС соседних частиц парный потенциал не является 
исчерпывающей характеристикой взаимодействия и применение ячеечного метода может 
существенно дополнить наши представления о свойствах ПКС. 
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